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We construct and test a quasi-perfect lattice action for staggered fermions. The construction 
starts from free fermions, where we suggest a new blocking scheme, which leads to excellent 
locality of the perfect action. An adequate truncation preserves a high quality of the free 
action. An Abelian gauge field is inserted in d = 2 by effectively tuning the couplings to a few 
short-ranged lattice paths, based on the behavior of topological zero modes. We simulate the 
Schwinger model with this action, applying a new variant of Hybrid Monte Carlo, which damps 
the computational overhead due to the non-standard couplings. We obtain a tiny "pion" mass 
down to very small (3, while the "77" mass follows very closely the prediction of asymptotic scaling. 
The observation that even short-ranged quasi-perfect actions can yield strong improvement is 
most relevant in view of QCD. 
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1 Introduction 



Most QCD simulations so far have been performed either using Wilson fermions or staggered 
fermions The latter formulation is especially useful in the chiral limit, because the remnant 
chiral symmetry U(l) <S> U(l) protects the zero fermion mass from renormalization. As a related 
virtue, its artifacts due to the lattice spacing a are only of 0(a 2 ), whereas they are of 0(a) for 
Wilson fermions interacting by gauge fields. 

It is now widely accepted that the above lattice actions should be improved, so that the 
lattice spacing artifacts are suppressed and coarser lattices can be used || . There are essentially 
two improvement strategies in the literature. In Symanzik's program Q the action is improved in 
orders of a. For QCD with Wilson fermions this has been realized on-shell to the first order on the 
classical level || and recently also on the quantum level || . Less work has been devoted to the 
improvement of staggered fermions, perhaps because the artifacts in the standard formulation are 
already smaller. However, S. Naik has applied Symanzik's program on-shell, where he improved 
the free staggered fermion by adding more couplings along the axes 0, and the Bielefeld group 
did the same by adding diagonal couplings |J. Furthermore, the MILC collaboration achieved 
a reduced pion mass by treating the gauge variable as a "fat link" |J. Finally, some work on 
improved operators has been done in this framework | |10| |, 

The other promising improvement scheme is non-perturbative in a and uses renormalization 
group concepts. It has been known for a long time that there are perfect lattice actions in pa- 
rameter space, i.e. actions without any cutoff artifacts [11]. More recently, it has been suggested 
to approximate them for asymptotically free theories as "classically perfect actions" |l2]|, which 



works very well in a number of two dimensional models [|Tj, 13, 14, |15|| , and it has also been 
applied in 4d pure Yang- Mills gauge theory |16|, [l?]] . 

For free or perturbatively interacting fields, perfect actions can be constructed analytically 
in momentum space. For Wilson type fermions this has been carried out to the first order in the 
gauge coupling in the Schwinger model [18] and in QCD [19|. A technique called "blocking from 
the continuum" was extremely useful for this purpose. One expresses all quantities in lattice 
units after the blocking, and sends the blocking factor to infinity. Hence the blocking process 
starts from a continuum theory, and it does not need to be iterated in order to identify a perfect 
action. 

For staggered fermions, a block variable renormalization group transformation (RGT), which 
does not mix the pseudo flavors and which does therefore preserve the important symmetries, 
has been suggested in Ref. [20]. It requires an odd blocking factor n. Iterating the n = 3 block 
variable RGT, a fixed point action, i.e. a perfect action at infinite correlation length, has been 



constructed [13, 21]. Also for staggered fermions blocking from the continuum is applicable |2^ 
This has been carried out for a general (flavor non-degenerate) mass term, revealing the intimate 
relation to the Dirac-Kahler fermion formulation in the continuum [23], and also including a 
suitable treatment of the gauge field [24]. Using the generalization to a flavor non-degenerate 
mass ]23[| , the spectral doublers inherent to the staggered fermion formulation might be treated 
as physical flavors in a QCD simulation. 

In the present paper, first the procedure of blocking staggered fermions from the continuum 
is revisited. We then discuss the optimization of locality, in the sense of an extremely fast 
exponential decay of the couplings in coordinate space. This property is crucial for practical 
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applications, because we have to truncate the couplings to a small number, which is tractable 
in simulations. Of course the truncation violates the perfectness, but for excellent locality this 
violation is not too harmful. 

We suggest a new blocking scheme, which we call "partial decimation" . It leads to a higher 
degree of locality than the usual block average method, i.e. to a faster exponential decay of 
the couplings in coordinate space. We then truncate the couplings to a short range by means 
of mixed periodic and anti-periodic boundary conditions, which are particularly adequate for 
staggered fermions. The excellent quality of the truncated perfect action for free staggered 
fermions is confirmed by spectral and thermodynamic considerations. 

We then proceed to the two flavor Schwinger model [25] as a test case. The Schwinger 
model is well-understood from continuum calculations, also in finite volume [^, [2?], |2^[. It 
shares important features with QCD, such as asymptotic freedom, confinement and topological 
quantum numbers with corresponding zero-modes of the Dirac operator. On the other hand, 
one expects other features of higher-dimensional gauge theories not to be well-represented by 
the Schwinger model, due to its super-renormalizability. 

The fermion gauge vertex is added "by hand" to the truncated perfect free fermion. We first 
insert the U(l) links between fermionic source and sink along certain shortest lattice paths. In 
a second step, we implement "fat links" for the paths consisting of just one link. The staple 
weight is determined effectively by minimizing the eigenvalues of the approximate topological 
zero-modes. For the pure gauge part, we use two actions, which are perfect resp. approximately 
perfect in d = 2. 

Of course, the ad hoc treatment of the vertex deviates from the systematic construction of a 
(classically) perfect action. However, in the staggered scheme any approximation or truncation 
of perfect action can cause errors of 0(a 2 ) at most, and the remnant axial symmetry protects 
the mass term from renormalization. By contrast, a strong mass renormalization is a severe 



problem for quasi-perfect fermions of the Wilson type p9, 30, 31 



Naively, the computational effort for simulations increases about linearly with the number 
of couplings in the action. In addition, physical effects, e.g. more exact zero-modes, may be 
a hurdle for the application of an improved action. We propose a way to damp this increase 
of computational costs, exploiting the freedom to design the Molecular Dynamic steps in the 
Hybrid Monte Carlo algorithm. We use a simplified action there, and the full quasi-perfect 
action in the Metropolis acceptance decision. We discuss the performance of this method for 
different values of /3, the crucial question being the acceptance rate. This method may also 
accelerate simulations with improved actions on parallel machines. It depends on a sufficiently 
good simplified version of the improved action, which is another plus for staggered fermions. 

For quasi-perfect staggered fermions, the dispersion relation and the scaling behavior of the 
mass spectrum are good down to (3 ^ 1.5. For these low values of /3, the masses for the n- 
and //-particle demonstrate a good scaling resp. asymptotic scaling behavior. In particular the 
scaling is even better than the one observed in Ref. [pi]], which uses truncated perfect Wilson 
type fermions, and 123 independent couplings parameterizing a classically perfect fermion-gauge 
vertex. Compared to that scheme, a drastically reduced number of couplings is needed here, 
and even more in the step to d = 4. Our results demonstrate that in fact a relatively modest 
number of couplings can yield a very powerful improvement. 

A synopsis of this work was presented in Ref. [33]. 
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2 The staggered blocking scheme 



The construction of perfect actions for free staggered fermions by blocking from the continu- 
um has been described in Refs. [23, ^3]. We start by briefly reviewing this procedure. The 



staggered blockspin fermions are defined in two steps. First we transform the Nf = 2 d l 2 fla- 
vors of continuum Dirac spinors ipa(x) (a: spinor index, b: flavor index) into the Dirac-Kahler 
(DK) representation given by (p(x,H). Q These functions originate from the representation of 
inhomogeneous differential forms 

5^ tp(x, H) dx H , dx H = dx^ 1 A... Adx^ , (2.1) 

H 

where H = {hi, ■ ■ ■ ,fJ>h}i /Xi < . . . < ^ is a multi-index. Transformation and inverse transfor- 
mation read 

ip(x,H) = * J^*^), 7* = W^-V*, (2-2) 

V N f ab 

= (2.3) 

^ f H 



Second, we introduce a coarse lattice of double unit spacing f = {y \ y^ = 2^}, which is 
a sublattice of V = {y \ y^ = n^}, with n^ra^ £ /Z. The fine lattice points y are uniquely 
decomposed as {fx is the unit vector in /^-direction) 

V = V + e H , e H = ^2 A • (2-4) 

tieH 

Thus the multi-index H{y) defines the position of a fine lattice point y with respect to the 
coarse lattice T. Now the blockspin variables ^(y) can be defined as averages of the component 
functions <p(x, H(y)), with a normalized weight U(x — y), JdxU(x—y) = 1, which is assumed to 
be even and peaked around x = y, 



= ll {y) * i dx U(x - y) ^(*) . (2.5) 

f ab 



This scheme has been proposed first in Ref. | 22| . Its peculiarity is that the staggered block 



centers depend on the multi-index H of the Dirac-Kahler component functions f(x,H). Block 
average (BA) means in this case average over the overlapping lattice hypercubes [y] = {x\ — 1 < 
(Xfj, — Vfj) < 1}. This scheme is given by n = Uba{x) = 2~ d for x G [y], and Hba(x) = 

otherwise. 

For the following calculation we diagonalize the lattice action using the staggered symme- 
tries. Here it is important that fine lattice shifts are no symmetry transformations. However, 
combination with site-dependent sign factors gives rise to the non-commuting flavor symmetry 
transformations |36|]. Therefore we replace ordinary Fourier transformation by harmonic anal- 
ysis with respect to flavor transformations and coarse lattice translations. We thus obtain a 

2 For the relation of the DK formulation of continuum fermions with s tag gered lattice fermions we refer to 
Ref. tea]. The relation to the block spin transformation is discussed in Ref. []23[. 
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modified momentum representation which intertwines Fourier transformation and the transition 
back from DK fermions to the Dirac basis, 

y 

*£(p) = E eiPy t£ W *(V) > P e fi = ] - tt/2, vr/2] d (2.6) 



(S is the Brillouin zone with respect to the coarse lattice). Inserting Eq. (2.5) we find 

<(p) = 4r n ^') <V) E ^ (p " p% 7^ fe) 7S y) * , (2.7) 

V ly f J n a >y y 

where tp^(p),H(p) denote the Fourier transforms of tp^(x),Il(x). The last sum can be re-written 

as 

E ^ E e iqeK t5 i% * = * d E % - *o E II ^ ^ * 

= A^E^-^S ^ • (2-8) 

We have used the orthogonality of the 7-matrix elements 

£tS ^V* = NfSaa'Sw , (2.9) 
H 

and H(l) is defined by H{1) = {fi | Lis odd}, H = H for /i even, H = {n \ (j, £ H} for /i odd. 
Finally Eq. ( |2.7D becomes (summation over double spin and flavor indices is understood) 

<{P) = y/Nf E n(p + 7 rZ) ^(p + ttQ , #(p + ttZ) = 7 fj° i#(p + *Q 7?^- (2-10) 

Note that the blockspin transformation is diagonal with respect to spin and flavor for continuum 
momenta within the first Brillouin zone B, yet not for all I ^ 0. 

We are now prepared to compute the perfect action for a RGT of the Gaussian type. Starting 
from a continuum action with a general mass term - which does not need to be flavor 
degenerate - the perfect lattice action S[&, Vl/] is defined as 



:»,»] = JvmjvfjVr, e Xp {-y^^(-g)(^ a ,g M + m 6 )^,(g)} 



x exp{ f * [ [¥*(-p) -^^n(p + ttO ^ a b (-p - TTi) ] r/»(p) 

+V b a(~p) [ <(P) ~ y/fy E n(p + 7rf) ^Cp + ^) ] 

+^("P)^(P)^(P)]} , (2-11) 
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where rj, r\ are auxiliary Grassmann fields denned on the same sites and with the same flavor 
structure as ^f. A non-zero term 



(2.12) 



"smears out" the blockspin transformation, as in Ref. [24|. This term is used to optimize locality 
of the resulting perfect action; it will be specified later on. The Gaussian integrals over ip,il> and 
77,77 can be evaluated by substitution of the classical fields, which leads to | 

dp 



S[V, *] 
with the lattice propagator 



Nf TT d 



B IV f 



m-p) G- x Zip) <>{p) 



(2.13) 



6(1) [~'i{-^) K laa'{P + ^l)^ + 'm d ]5 ddl #(i)t 



(p + 7rZ) 2 + m\ 



7 d T • (2-14) 



Note that G is flavor diagonal, G bb ' = G b 5 bb ' , because 7^)7^ is diagonal iff 7^ is. In 

particular, for a degenerate mass term the adjungation with 7^) is trivial, and the lattice 
propagator is proportional to in flavor space. We define 



G h m >(p) 



(2.15) 



hence Q 1 ^ , Qq become 



2 ' 



Q b ( P ) = D b (p) + j2^(p + ^) 2 ^J2J2^( l )^i^* 



my 



N 



f y K 



(p + 7t/) 2 + m?, 



(2.16) 
(2.17) 



In case of a non-degenerate mass mj, the sums J2y J2k cannot be contracted according to 
Eq. (^9|), due to the sign factor ex(l) = Y[u&k(~ ^) lv ■ However, in the degenerate case m& = m, 
with flavor independent smearing terms D b = D, we simply obtain 



Q» = D^{p) + J2 n (P + ^) : 



(_l)^(p + 7r /) M 
(j) + irl) 2 + m? 



Qo(p) = D ( P ) + £ it> + ^) 5 



(p + 7rZ) 2 + m 



2 ■ 



(2.18) 
(2.19) 



The perfect action in real space arises from Eq. ( 2.13 ) inserting the momentum representation 
given in Eq. (2.6). After some 7-matrix algebra (2*|], we arrive at 



= ^^(y)m(y,y')^(2/ / ) , 

v,v' 

m(y,y') = ^2p K (y')p(y -y',y') M K (y -y') . 



(2.20) 
(2.21) 



K 



We ignore constant factors in the partition function. 
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Corresponding to a lattice propagator diagonal in flavor space, the sum over K runs over multi- 
indices K £ £> with diagonal j K , in the Weyl basis T> = {0, 12, 34, 1 2 3 4}. The sign factors 
p(z, y) = p(H(z),H(y)) arise from j h j k = p{H, K) j HAK , where HAK = (H U n if), 

and p K (y) is given by p(H(y), K)p(K, H(y)). By symmetry, the only non-zero contributions to 
M K (y) are 

M*(y) = ip(p,K) [ M«(p) for H(y) = p,AK , (2.22) 

M <f (y) = / * e ~ w M o x b) for # G/) = K . ( 2 - 23 ) 

w Ml m5, w - ^ ?7 - S! ^| f |L s _. „ 24) 

The flavor degenerate case leads to vanishing components for K ^ 0, and (with M^ = M^$) 
we simply obtain 



M,,o(p) = ^ n C l^ n , V2 • (2.25) 



It has been proven in Ref. |2J| that the couplings given by the fermion matrix m(y, y') are 
local, i.e. they decay faster than any power of \y — y'\. For that, certain periodicity properties 
apply, which translate into (for simplicity of notion, p denotes either p or 0, and KAO = K) 

Qlip) = Y.-&Qt(P)< <£&> + «*> = I nil^ til Ifr, ■ < 2 ^> 



Kev 



Qf (p) for v £ KAp 



Again, we sum over diagonal 7-matrices only. It is provided that the corresponding requirements 
are met for the smearing terms Di,-, within Q~(p), see below. In consequence, the fermion 

matrix components Mjf{p) obey periodicity conditions analogous to Qf[(p), and the integrands 
of Eqs. ( [2.22 , 2.23j ) are periodic with respect to the Brillouin zone B and analytic in a strip 



around the real axis. This implies locality of the perfect action. 

The coupling of even and odd lattice points is due to the Mjf components of the fermion 
matrix; the M$ components couple even-even and odd-odd. We add without proof that even- 
odd decoupling of the Hermitian matrix m can be shown in any even dimension d with 
arbitrary (non-degenerate) mass terms for truncated versions of the perfect fermion matrix 
m(y,y'), see Ref. |23| for d = 2. This is a useful property in simulations with Hybrid Monte 
Carlo algorithms. However, with (non-perfect) coupling to a gauge field and non-zero mass term 
(i.e. with even-even, odd-odd as well as even-odd couplings), this is not true in general. |^ 

3 Optimization of locality 

In the blocking scheme described so far, there is quite some freedom left. In particular, we 
may use averaging functions different from Hba, and we can choose the smearing term D in 
Eq. (Imp . In both cases we aim at optimization of the locality in the resulting perfect action. 



4 Since even-odd decoupling of m is a perfect property, it could be imposed as a construction requirement 
for an approximately perfect fermion-gauge vertex. For m = 0, even-odd decoupling of m is guaranteed (m 
only couples even with odd sites) . 
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Let us first discuss the averaging scheme. In Eq. (2.11) we implicitly assumed the same 
blocking of t/j and t/>, given by the weight function resp. its Fourier transform H(p). Now 
we consider the case of a block average II = TLba for ip only, while i/S (ip) is put on the lattice 



by decimation, II(x) = 5(x). Thus we obtain a sing le factor II (p + irl) in Eq. (Og) with 



n(p) = ii BA { P ) = n — ' P» = sin ^ • (3-1) 



In case of a 5 function RGT (D = 0), this means to identify the averaged continuum and lattice 
2-point functions 

( <Kv) W) > = I dx ( <P(?,H{y)) <f(y',H(y')) ) . (3.2) 

■%] 

Due to translation invariance it does not matter whether we average source or sink of the 
continuum expression, or whether we allocate the space directions to be integrated over to 
source and sink in some way. (The last point of view may be used to make a closer contact to 



the construction of staggered fermions from DK fermions in the continuum [35|, as discussed 
in Ref. ]23|] .) We call this blocking scheme partial decimation. For the 2-point functions every 
space direction is integrated over once; therefore we do not run into the difficulties arising for 
blockspin transformations with complete decimation, which do not have a corresponding perfect 
action. 

For both blocking schemes, block average for ip and ip (BA) and partial decimation (PD), we 
now want to optimize locality of the couplings by making use of the smearing terms in Eqs. ( 2.16| , 



2.171 ). As an optimization criterion it has been suggested to require that in the effectively Id 



case - with momenta p = (p\, 0, . . . , 0) - the couplings are restricted to nearest neighbors as in 
the standard action [19, ^]J^| In the degenerate case we require 



*7iMi(pi J 0,... ) 0) + Mo(pi,0,...,0) = /(m)[ipi7i + m], (3.3) 
with m\ m= Q = and /(0) = 1. Our ansatz for the Gaussian smearing term reads 

D n(p) = c{m)p^ D (p) = a(m). (3.4) 



Requirement (|3.3|) can be fulfilled in both blocking schemes we are considering, if we specify the 
RGT as follows 

CBA(m) = [coshm — 1] / (2m) 2 , as^(m) = [sinh(2m) — 2m] j (2m) 2 , 

cpu(m) = 0, ap^(m) = [cosh m — 1] / (2m) . (3-5) 

In both cases, we obtain m = sinhm. For m = a non- vanishing static smearing term a(0) 
would explicitly break the remnant chiral symmetry in the fixed point action. Therefore, the 
static term should vanish for optimized locality, as it does in both cases. Furthermore, in the PD 
scheme the chiral limit is optimized for locality by a simple 5 function RGT. As an advantage of 



It has been shown in Ref. |S7j that this criterion does optimize locality in d = 4 for scalar fields over a wide 
range of masses. 
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this property - which is not provided by the BA scheme - there is a direct relation between the 
n-point functions in the continuum and on the lattice. In addition, the extension to interacting 
theories might involve numerical RGT steps in the classical limit, which also simplify in the 
absence of a Gaussian smearing term. Finally, a non-vanishing term c(m) causes complications 
if one wants to include a gauge interaction in the RGT, but the PD scheme avoids such problems. 

The decay of the couplings m(x, 0) = m(x) in the massless case for d = 2 and d = 4 is shown 
in Figures 1 and 2. We see that the PD blocking scheme works better. In the non-degenerate 
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Figure 1: The decay of couplings in (1,0) direction for d = 2 in the massless case: BA without 
smearing (diamonds), BA optimized (squares), PD (triangles). 

case, the simplest ansatz for a smearing term is (corresponding to the staggered fermion action 
with non-degenerate mass |38|]) 

Dfo) = cp^ , D b (p) = Ibh II cos ^ • (3-6) 

KeV ntK 

In this case, we obtained - by means of numerical optimization - a similar decay of couplings 
as with degenerate masses, yet no strict 1-dimensional ultralocality. For d = 2 it appears that 
a non-degenerate parameter au does not improve the coupling decay significantly, so we 
worked with a\2 = 0. Again the PD scheme, where we assumed cpd = 0, leads to a more local 
action, see Figure 3. 

4 Truncation effects 

The litmus test of any perfect action is given by the truncation effects in a practicable number 
of remaining couplings. These effects are minimized by maximal locality. Yet, the truncation 
scheme itself may have some impact on the truncation errors. An elegant procedure has been 
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Figure 2: The decay of 4d couplings in various directions in the massless case: BA with optimal 
smearing term and PD. We see that the latter couplings decay faster. 



proposed in Ref. [gjj for Wilson fermions. First the perfect action is constructed on a small lattice 
volume N d by restricting the momentum components to the discrete values = 2nn ll /N € 
] — 7T, 7t] 3 n M € Z,. Typically, one chooses N = 3, and the resulting couplings are then used 
in a large volume too, where they are not exactly perfect any more. This truncation scheme 
has the virtues of automatically correct normalizations, and a simplification of the numerical 
evaluation as opposed to a truncation in coordinate space. Furthermore, the mapping to the 
corresponding truncated perfect action in a lower dimension remains exact; for instance, one can 
reproduce the effectively Id nearest neighbor action starting from d > 1 by summing over the 
extra dimensions, which provides a sensitive test of the numerical accuracy. Provided a good 
locality, this method works well for Wilson-type fermions, gauge fields |2{| and scalars [p7|1 . 

For staggered fermions, only a crude truncation in coordinate space has been applied so 
far One includes the coupling distances ±1, ±3 in the [i direction, and 2,0,-2 in the 

non-/i directions. This involves more couplings than the N = 3 truncated Wilson-type fermion 
(called "hypercube fermion"), whereas the number of degrees of freedom is the same in both 
cases. However, although significant improvement has been achieved compared to the standard 
staggered formulation, the quality of its spectral and thermodynamic properties did not reach 
the level of the "hypercube fermion". In order to arrive at results of a similar quality, we 
adapt the above truncation procedure - together with the PD scheme - to the case of staggered 
fermions. 

We treat the components of the fermion matrix M~f(p),jl = /U,0 separately, as they show 
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Figure 3: The decay of couplings in (1,0) direction for d = 2 with masses mi = 0, ni2 = 0.5: 
BA without smearing (diamonds), BA optimized (squares), PD (triangles). 



different behavior under reflections, 

*,K, ^ f ~ M ti (?) for^GifA/i 

M, ( Pl . . -ft,,ft, +1 ... Pd) = | for ^ . (4.1) 

Again we use the short-hand notation KlS.fi = K for /x = 0. Truncation is achieved by discrete 
Fourier transformation and a discrete support given by cjv(y^) = 1,1/2,0 for |y M | <,=,> N, 
respectively, 

= n^f^ x: e ~ ipy Mjx(p) . (4.2) 

It is the set of discrete momenta 6~^y which depends on the reflection properties given by \x, K. 
We choose mixed periodic and antiperiodic boundary conditions, 

nK I (n + l/2)n/N for v G i\A/i , ,„ /ri , ,. „, 

P e ^ = ( ^ for „ ^4 , n G ft, € ] - vr/2,vr/2] . (4.3) 

It is easily verified that the transformation of the truncated components back to momentum 
space reproduces the perfect values at p G B^. N , 

MgM = £ eiPV M *n(v) = M f (P) fOT ^B^. (4.4) 

The components M~f(y) inherit the reflection behavior of M~" (p) described in Eq. ( f4.1|) . There- 
fore, treating them as periodic functions in y-space by discrete momenta p u = nn/N,n G 2Z 
would make them vanish artificially on the boundaries y u = ±N for v G KAjl. We avoid this 
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0.5 1 1.5 

Pi 

Figure 4: The massless spectrum in d = 2. We compare standard staggered fermions (crosses) 
and perfect spectrum (full line) with perfect actions truncated with N = 3 and mixed periodic 
boundary conditions: block average (dotted), partial decimation (dashed). The last case is also 
plotted for truncation with periodic boundary conditions (dashed-dotted) . 



defect by the above choice of discrete momenta. It pays off by a drastic reduction of truncation 
effects, see Figure || for the 2d spectrum in the massless case using partial decimation. The 
values of the perfect couplings, truncated by mixed periodic boundary conditions (as described 
above) for N = 3, are given in Table |]. 

Furthermore, the truncation effects are significantly stronger for even truncation distances 
N; in particular, for N = 2 and N = 4 the energies become complex-valued at large momenta. 

For N = 3 - which is still tractable in numerical simulations - the amount of improvement 
is already striking. We compare the spectra for d = 2 in the PD and BA blocking scheme 
in the massless case (Figure ||) and for non-degenerate masses mi = 0, ni2 = 0.5 (Figure ||). 
The couplings derived from the PD blocking scheme turn out to be better, in agreement with 
the higher degree of locality observed in the (untruncated) perfect action. We see that the 
improvement is still good in the case of non-degenerate masses. In this case, we optimized the 
smearing parameters numerically, as pointed out in Section ||[ The values are clba = 0.120, 
cba = 0.125, apD = 0.082, cpd = 0. The standard staggered fermion results for non-degenerate 
masses are calculated with the action proposed in Ref. |3Sj| . We emphasize that the spectra of 
(untruncated) perfect actions are indeed perfect, i.e. identical to the continuum spectra (up to 
periodicity, which is inevitable on the lattice). Hence the spectrum reveals directly the artifacts 
due to truncation. 

Figure 6 compares the massless spectra in d = 4, and we see that the qualitative behavior 
observed in d = 2 persists. We compare the standard staggered fermion, and the optimized BA 
and PD fixed point fermions, both truncated by N = 3 mixed periodic boundary conditions. 
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Figure 5: The spectrum in d = 2 for m\ = 0, 771,2 = 0.5. Standard staggered (crosses) and 
perfect spectrum (full line) are compared with perfect actions truncated with N = 3 mixed 
periodic boundary conditions: block average (dotted), partial decimation (dashed). 



We see again that the PD scheme is superior. Furthermore, we show for comparison also 
the dispersion relation of a Symanzik improved action called "p6" from Ref. ||. Symanzik 
improvement by additional couplings along the axes (Naik fermion, 0]) only yields a moderate 
quality [24], but the Bielefeld group suggests a number of actions, where Symanzik improvement 
is achieved by diagonal couplings. The p6 action is the best variant among them, and also Figure 
6 confirms its excellent level of improvement. However, it is not obvious how that formulation 
can include a general mass term in a subtle way. 

As a further test in d = 4, we consider two thermodynamic scaling quantities. The first 
is the ratio P/T 4 (P: pressure, T: temperature). According to the Stefan-Boltzmann law, 
this ratio is 77r 2 /180 for massless fermions in the continuum, and a lattice action with many 
discrete points iVj in the temporal direction will asymptotically reproduce that value. However, 
the speed of convergence, and in particular the behavior at small N t , depend on the quality of 
the action. In contrast to the spectrum, this ratio is not even exact for the fixed point action 
in Eqs. ( 2.18 , 2.19 ), because of the "constant factors" that we ignored when performing the 
functional integral in the RGT Eq. ( |2.1rj ). Such factors may depend on the temperature, so 
with respect to thermodynamics our action is not fully renormalized |29|] . However, it turns out 
that the unknown factor is very close to 1, except for the regime of small Nt (about Nt < 6), 
which corresponds to very high temperature. So the main issue is again the contamination due 
to the truncation. 

In Figure 7 we compare this thermodynamic scaling for a variety of staggered fermion ac- 
tions at 777 = 0, and again the PD scheme turns out to be very successful. A similar level of 
improvement can be observed for the p6 action, which is also here by far better than the Naik 



13 



d = 4 






d = 2 








PD 


BA 




PD 


BA 


(1,0,0,0) 


0.348194 


0.331558 


(1,0) 


0.4391168527 


0.433150 


(1,2,0,0) 


0.020490 


0.022963 


(1,2) 


0.0304416413 


0.033425 


(1,2,2,0) 


0.002240 


0.002435 








(1,2,2,2) 


0.000247 


0.000180 








(3,0,0,0) 


0.007609 


0.011721 


(3,0) 


0.0052123698 


0.008022 


(3,2,0,0) 


-0.000216 


-0.000319 


(3,2) 


-0.0026062073 


-0.004011 


(3,2,2,0) 


-0.000384 


-0.000605 








(3,2,2,2) 


-0.000214 


-0.000317 









Table 1: The couplings of the perfect action for massless staggered fermions, constructed from 
the "partial decimation" (most successful) and from the optimized "block average" scheme, and 
truncated by mixed periodic boundary conditions. The couplings are odd in the odd component, 
and even in all other components. Among the latter there is also permutation symmetry. The 
2d PD couplings (given to higher accuracy) are used in the Schwinger model simulation reported 
below. 

fermion. We conclude that a good improvement method should in any case include diagonal 
couplings, which are far more promising than additional couplings on the axes (just consider 
rotational invariance, for example). Thus we can expand the range, in which a practically ac- 
curate continuum behavior is observed, by about a factor of 4 compared to standard staggered 
fermions (see Figures 4, 5, 6 and 7). 

As a second example for thermodynamic scaling of free massless fermions, we set T = but 
switch on a chemical potential [i. The ratio of the "baryon number density" tlb divided by /x 3 
is our second scaling quantity. If fi approaches 0, this ratio converges to its continuum value 
2/(97r 2 ) for any lattice action. For increasing \i the lattice artifacts become visible, and we see 
again that the truncated perfect staggered fermions, especially those based on the PD scheme, 
stay much longer in the vicinity of the continuum value, see Figure 8. Again the practically 
accurate regime is extended by a factor ^ 4. 

5 The fermion- gauge vertex function 

We construct a quasi-perfect fermion-gauge vertex using the truncated perfect couplings of the 
free fermions based on the PD scheme, given in Table 1. Without a mass term, only even-odd 
couplings are present. The coupling to the gauge field is achieved by the insertion of U(\) 
parallel transporters on shortest lattice paths between the sites of ^ and Where several 
shortest ways exist, we use only staircase- like paths, as illustrated in Figure ||. For and on 
nearest neighbor sites we also allow for the shortest detour (staple) with some weight factor w, 
leading to "fat links" . Therefore, the fermionic part of the action takes the form 

Sf[V,V,U] = Y.Hp^ x ) { °i.o (l-2w) *(x) U^x) *(x + e M ) 

X fl 
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Figure 6: The dispersion relation for various types of 4d massless staggered fermions along 
(1,1,0). We compare the standard action, a Symanzik improved action called "p6" from Ref. 
[8] and the optimized BA and PD hxed point actions, truncated with mixed periodic boundary 
conditions. 

+ c lfi wJ2^(. x ) U v {x)U^(x + e v )U- v (x + e„ + e M ) ^(x + e^) 
+ ci )2 U u (x)U^(x + e u )U u (x + e u + e M ) V(x + + 2e„) 

+ c 3 ,o * (a) U^U^x + e^U^x + 2e M ) (x + 3e M ) 
+ c 3 ,2 51 *0*0 U^{x)U v {x + e^L/^x + + e^fT^x + 2e^ + e„) 

x [^(x + 2e M + 2e u ) V(x + 3e^ + 2e„) } , (5.1) 

where we sum over p = ±1, ±2, z/ = ±/2, and fl is the direction perpendicular to p. For positive 
p the sign factor reads 

p{n,x) = (-If, n = , (5.2) 

and negative p follow from = — e^, p(—p,x) = —p(fi,x), f7_ M (x) = J7 M (x — e^)*. 

To optimize the fat link weight w, we study the behavior of the fermion matrix in a gauge 
field background of topological charge Qt op = 1- In the continuum limit this matrix should 
have two zero eigenvalues (corresponding to the two flavors), and we choose w such that we 
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Figure 7: The Stefan-Boltzmann law for standard staggered fermions and for the "p6 action", 
compared to the truncated perfect fermions for the optimized BA and PD scheme. 



reproduce this property on the lattice as accurately as possible. To this end, we generate a 
sample of quenched configurations to a fixed topological charge and a fixed (large) value of 0. 
In these configurations there are two fermion matrix eigenvalues of particularly small absolute 
value. We tune w in order to minimize their average, Ao- This yields w = 0.238. If we fix 
the physical charge and size as e c , L c , so that = l/(ae c ) 2 , L = L c /a, then this choice of w 
leads in the continuum limit to an improved behavior of Ao oc a 3 , in contrast to Ao oc a 2 for 
w = 0, see Figure [l(]. Of course, the determination of w is only approximative, because of the 
use of quenched configurations. The corresponding unquenched study should be also feasible. 
However, in physical theories like QCD, it is much more expensive to tune the parameters with 
unquenched configurations. 

6 The perfect pure gauge action 

We set the pure gauge part of the action either to the standard Wilson plaquette action SW[^1 
- which is perfect with respect to a somewhat complicated RGT [l^] - or to an approximation 
of another perfect action, which is of the Villain type. The construction of that action is closer 
to the PD scheme at m = 0, and its approximation is given by 

S V [U] = /?£F(x) 2 , e iF ^ = P(x), F(x)E(-7r,7r], (6.1) 

x 
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Figure 8: The scaling ratio tibIp? at T 
function of the chemical potential \i. 



for different types of free staggered fermions, as a 



P(x) = Ui{x)U 2 (x + ei)Ui(x + e 2 )* U 2 (x)* . (6.2) 

The action SV[£7] can be constructed in d = 2 by requiring the identity of correlation 
functions of gauge invariant quantities in the continuum and on the lattice. This parallels the 
treatment of the fermions that we are using here (we recall that in the PD scheme we block 
massless fermions by a 5 function RGT). We work with the (real) phases Au(y) (non-compact 
gauge fields) of the parallel transporters U^(y), and define a plaquette phase as 

F(y) = + A 2 (y + ei ) - A^y + e 2 ) - A 2 {y) . (6.3) 



It is equal to the lattice field strength F{x) from Eq. (|6.lD modulo 2n. In the continuum we 
define the phase of parallel transporters around plaquettes [y] corresponding to the above lattice 
quantity by 



f(y) = / a^dx^ , (6.4) 
Jd\y] 

where is the continuum gauge field. We require the equality of the continuum and lattice 
two-point functions 

lattice — continuum • 

(6.5) 

This leads to the definition of a perfect action 



e S[A] 



J Va,\[5{F(y) - f(y))e^I dx ^ 2 , (6.6) 
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Figure 9: The parallel transporters representing the 5 independent couplings. 



where f(x) is the continuum field strength. We assumed the perfect action to depend on the 
variables F y only, i.e. on all closed loops except Polyakov loops, like the continuum pure gauge 
action. 

For the continuum action we may decompose the field strength f{x) in averages f(y) over 
the plaquettes [y] and fluctuations fh(x) (with f, idx fh(x) = 0) 



Jdxf(x) 2 =J2j dx(f(y)+f h (x)) 2 =J2Ky) 2 + Jdxf h (xf 



(6.7) 



Furthermore, we decompose the continuum configuration integral into an integral over all gauges, 
over the torons (the constant part of the gauge field a M (x) , see Section ||) , an integral over the 
field strength configurations [/(x)], and an integral over the vortex configuration [n y ] for each 
plaquette, 

n y = (f(y) - f(y))/(2n) . (6.8) 

The deviations from Stokes' theorem arise for topologically non-trivial gauge fields. For the 
[/(l)-topology on the 2d torus T2 it would be sufficient to allow for vortices at one specific point 
x. However, for a local formulation, we may decompose T2 into patches corresponding to all 
plaquettes [y], see Ref. |39[| . 

The integral T>[f(x)] factorizes with respect to f(y) and fh(x), and we are left with 



-S[A] 



EE / d f(y) - f(y) - 2 ™y) e 

y n v 



-Pf{yf 



(6.9) 



which leads to 



S[A] = $>(^) , s{F y 



In 



-P(F(y)-2Trn y ) 2 



my? 



(6.10) 



Remember that F(y) is the projection (mod 2tt) of F(y) into the standard interval (— vr, ir]. The 
above approximation, Sy in Eq. (6.1), neglects only terms of order e _47r ^, if F y is not near an 
exceptional value ±ir (which is strongly suppressed at > 1). 
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Figure 10: The average lowest absolute eigenvalues Xo(P) for w = (+) and w = 0.238 (x ). 
Physical size and coupling are fixed so that 1? / '(3 = 36. The fits correspond to Ao oc 1//3 for 
u> = 0, and A oc 1//3 3 / 2 for ty = 0.238. 



Note that the above block scheme relies on the fact that - due to Stokes' theorem - blocking 
the 1-form a on links is consistent with blocking the 2-form f = da on plaquettes modulo vortices, 
see Eq. (|6.8j). In this sense, we proceed in analogy to the derivation of staggered fermions from 
DK fermions [35], see also the discussion in Ref. [23]. 

As a last remark on the pure gauge part, note that it is not straightforward to follow the 
above procedure for U(l) gauge theories in d > 2. It is easy to derive the condition for the 
lattice two-point functions corresponding to Eq. (|6.5|), but the lattice action cannot be derived 
by a direct inversion of the two-point function, since the plaquette variables are not independent 
any more. 



7 An adequate Hybrid Monte Carlo algorithm 

The feasibility of simulations with a large number of couplings in the action is a major obstacle 
against the use of quasi-perfect actions. The computational cost of such calculations is mainly 
given by frequent solutions of linear equations with the fermion matrix, which in turn need a 
large number of matrix multiplications. If the fermion matrix m^]^] is at hand, one matrix 
multiplication means n coup i x V complex multiplications, where n coup i is the number of path 
couplings to a fixed spinor ^(x) (or \P(a;)), and V is the lattice volume. For our truncation, one 
has n coup i = 24 in d = 2, compared to n coup i = 4 for massless standard staggered fermions, see 



19 



Eq. (5.1). In spite of the relatively modest overhead, saving the fermion matrix (n coup i x V com- 
plex numbers) - before any solution of a linear equation - requires a large storage space, which 
can be a serious problem in d = 4. On the other hand, any complicated vertex structure, i.e. 
additional parallel transporters between ^(x) and *$>(y) on fixed sites, do not require additional 
computational costs in leading order, since the structure of "hyperlinks" remains unaltered, see 
also Ref. p^ ]. An example for such additional parallel transporters are the staple terms, which 
we used to build the specific hyperlink called "fat link" . 

To suppress the increase of computer time needed, which in our case would naively amount 
to a factor of 6, we make use of the structure of the Hybrid Monte Carlo algorithm |Ifl]. It 
consists of a proposal, which is derived by a numerical integration of Hamiltonian equations in 
an artificial time (the Molecular Dynamics (MD) step), and a Metropolis acceptance decision. 
Detailed balance is guaranteed by a proposal symmetry within the MD steps, and the Metropolis 
decision with respect to the complete action S. The role of the MD proposal is to push the system 
far in configuration space without deviating much from the hyperplane S = constant. Thus the 
acceptance probability P a in the Metropolis decision can be kept at P a ~ 0.5 for large integration 
intervals in the artificial time, and therefore large moves of the system in configuration space. 
However, we are free to choose another action S' to be held approximately constant. This will 
lead to smaller acceptance probabilities, therefore the difference of S and S' should better not 
be too large. 

We choose the standard action in the MD steps, including the fat link parallel transporters 
for w 7^ 0. The gauge action is not altered. Figure 11 shows the behavior of the acceptance 
probability P a ((3)'- it decreases as one approaches the strong coupling regime. This may be 
expected, since for lower 0, S and S' really describe different physics, namely approximate 
continuum behavior for the perfect action S and strong coupling behavior for S' . With the fat 
link included, the strong coupling region is pushed down to lower values of (3. Therefore we 
obtain useful probabilities of P a also for lower /3-values, down to (3 « 1.5. 

The relation of computer time needed for the MD proposal and for the Metropolis decision 
is approximately given by the number of MD steps needed to achieve a sufficiently precise 
integration of the Hamiltonian equations in the artificial time. In our case we use a random 
choice between 6 and 12 MD steps corresponding to an artificial time interval of 0.25 to 0.5. 
Since the full fermion matrix multiplication needs a factor of 6 more floating point operations, 
this leads to an overall factor for the computer time of « 1.5 with the above method, compared 
to 6 for MD proposals using the full quasi-perfect action. This must be compared with the 
additional autocorrelation time due to the decreased P a , see Figure 11, which shows that there 
is still a substantial gain left. We find it worthwhile to check, whether it is possible to proceed 
in a similar way in d = 4 models. This will depend on the much stronger increase in the number 
of couplings for quasi-perfect actions, on the number of MD steps in these models, and on the 
behavior of P a . 

There are other effects, which can increase the autocorrelation time when one approaches the 
continuum physics with improved actions. Particularly in the Schwinger model, the approximate 
zero-modes will become sharper, and for any observable with zero-mode contributions a special 
treatment becomes necessary. For instance, one may measure the contributions from different 
topological sectors separately, and determine their relative contribution after a re-weighting of 
the sectors in question (in the sense of a discrete multi-canonical simulation |4l]]), as outlined 
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Figure 11: The acceptance rate P a {0) on a 16 x 16 lattice. The actions involved are: standard 
staggered fermions with gauge action S\y in the MD step and free perfect in the Metropolis step 
fe); standard staggered with w = 0.238 and gauge action Sy in the MD step and quasi-perfect 
in the Metropolis step (x ); standard staggered fermions and Sw in both cases (+). (The terms 
'free perfect' and 'quasi-perfect' are defined in Section 8.) 



in Ref. |^2|]. Alternatively, one could use a Polynomial Hybrid Monte Carlo algorithm f43f] , 
recently proposed to deal with zero-mode contributions. In both cases, the problem of effectively 
changing the topological sector remains. For our case, we used instanton hits matched to the 
present topological zero-modes, see Ref. [fl4|]. 

We want to mention another two numerical problems. First, the calculation of connected 
contributions (see Section 9) requires the calculation of propagators (m^ 1 ) with sources at all 
lattice points x. Doing this for each configuration would need a multiple of the time needed to 
generate these configurations. One often uses methods like the noisy estimator to circumvent 
this problem. In our case, however, we decided to really calculate all propagators, but not for all 
configurations, since there are large autocorrelation times in the fermionic observables (probably 
due to topological quantities with large autocorrelation times). This makes it also possible to 
average over the propagator source for disconnected contributions, and thereby decrease the 
statistical fluctuations of the single measurements. 

Secondly, it is difficult to evaluate small masses from correlations with large relative errors. 
For artificial data of a decay c(t) = cocosh(m(t — L/2)), modified with uncorrelated Gaussian 
fluctuations 5c(t)/c(t) = e, m < e, we found our fitting procedure to yield results for m in 
a single cosh- fit, which were not resolved from zero despite of error bars of order 1CP 7 . We 
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therefore trust in decay masses only, if they are resolved from zero by their error bars. This is 
a problem for the evaluation of the 7r-mass, where similar effects appear, in particular for small 
fit intervals. 

Altogether we generated ~ 3 ■ 10 5 configurations with maximal autocorrelation times up to 
100 for topological observables, and up to 25 for other gauge field observables. For fermionic 
observables we measured after steps of 100 configurations. Even after that reduction we found 
autocorrelation times up to 3. 



8 Results for the torons 

Here we present the results for the free truncated perfect action for staggered fermions 

• with the minimal vertex (w = 0) ("/ree perfect'), together with the standard pure gauge 
plaquette action Syy, 

• and for our extended approximation of a perfect action (quasi-perfect: wfl = 0.238, gauge 
action Sy). 

First, we investigate the behavior of the toron part of the gauge field, i.e. its constant mode 
on the 2d continuum torus f\ 



tn = -j—r- / dx aJx) . (8.1) 

It is convenient to define a dimensionless quantity u M = (eL M /2-7r) t^. Under large gauge trans- 
formations it transforms as — > + n M , £ /Z. Its effective fermion-induced action F(u) 
(the I l°g °f the fermion determinant) decouples from all other parts of the gauge field, and is 
independent of the gauge coupling e, 

T(u) = 2 In 1 0i (u 2 + irui) \ - 2-ktu{, t = L 2 /L x , (8.2) 

where 0i is the Jacobian 0-function. For details of the calculation in the two-flavor case we refer 
to Ref. |27]. 

We therefore expect an approximately perfect behavior of toron-dependent observables al- 
ready with the free perfect fermion. We measured the Fourier coefficients c nijn2 of the toron 
distribution 

c nii1l2 = J d 2 u e2 r (") cos(27ruirai) cos(27rn2n2) , (8-3) 

and results for various (3 and lattice sizes are given in Table H There is no dramatic deviation 
from the continuum results within small errors, even with the vertex based on shortest paths 
only. This is clearly in contrast to the standard staggered fermions. 

6 For a lattice version of the average in Eq. (|g.l|), the 27r-ambiguity of the phases A ti (x) of the link variables 
U^{x) have to be resolved by fixing the gauge invariant lattice field strength F(x) €_(— 7T, 7r], see Eq. ( |6.l[ ), and 
by the requirement Fix) — A\(x) + A2(x + ei) — Ai(x + e^) — Ai(x), see Refs. WA ft5[. 
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6xl6-lattice, /3 = 10 


Cl,0 


C2,0 


co,i 


C0,2 


standard 
free perfect 
quasi-perfect 
continuum 


-0.7619(9) 
-0.7441(7) 

-0.7435 


0.3415(16) 
0.3084(15) 

0.3079 


-0.0061(15) 
-0.0079(11) 

-0.0075 


0.0014(11) 
0.0009(8) 

0.0002 


16xl6-lattice, = 3 


Cl,0 


C2,0 


CQ,1 


C0,2 


standard 
free perfect 
quasi-perfect 
continuum 


-0.299(8) 
-0.301(10) 
-0.322(9) 
-0.322 


0.047(5) 
0.038(8) 
0.052(6) 
0.043 


-0.290(7) 
-0.312(10) 
-0.322(9) 
-0.322 


0.040(5) 
0.050(6) 
0.038(7) 
0.043 



Table 2: Some toron expectation values (defined in Eq. (8.3)) in the standard, free perfect, and 
quasi-perfect case, compared with the continuum results. 



9 The "meson" spectrum 
9.1 Observables 

Finally, the "meson" spectrum is derived from the correlation functions of point-like and one-link 
lattice bilinears corresponding to the isovector and isoscalar current (see Ref . |42|] ) . 



isovector: 




-(*) 


isovector: 




-(*) 


isoscalar: 




-(*) 



: ("If 2 V(X)V(X) , 

: [f + ei ) _ + ei )*( X )] , 

: #(x)*(x + ei) + #(x + ei)*(x)] . (9.1) 

The isospin / = ±1/2 labels the two flavors of Dirac fermions corresponding to staggered 
fermions for d = 2. From the bilinears M a (x) ai(T2 , a = 0, 1—, 1+, we derive correlation functions 
at fixed spatial momentum pi in the usual manner, 



c. 



, Tia2 (t; Pl ) = £e**«i ( M« ia2 ( Xl ,t) M« 1(T2 (0,0) ) . (9.2) 



Here, it is helpful to consider the mechanisms which determine the current correlation functions 
in the continuum. There are three possible contributions to fermionic two-point functions in the 
2-flavor Schwinger model [ 26| , 27]. The connected part (with propagators from x to y and back) 



is the only non-vanishing contribution for the isovector current correlation. It decouples from 
the gauge field, except from its constant (toron) part. The corresponding intermediate particle 
is massless (we call it "pion"). For the isoscalar current correlation, the disconnected part 
(propagators from x to x and from y to y) also contributes. It can be treated by gauge invariant 
point splitting, which introduces an (anomalous) coupling to the gauge field. Combined with 
the disconnected part, the correlation is then determined by an intermediate massive particle 
(the "ry-particle" ) . 

The current correlation functions receive only contributions from gauge fields with zero 
topological charge. Gauge configurations in non-trivial sectors are suppressed by topological 
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zero-modes, and the currents do not couple to the poles in the propagators which arise from 
those zero eigenvalues. In contrast, the correlation functions of scalar and pseudoscalar densities 
pick up zero-mode contributions, which arise from the cancellation of these poles with the zeros 
of the fermion determinant. For a discussion on how these mechanisms are realized on the lattice 
with staggered fermions we refer to Ref. p2| . 

In a given time slice we cannot fix the quantum number Thus the corresponding corre- 
lation functions obtain contributions from a 2 = ±, i.e. with and without a sign factor ( — 1) X2 . 
In this way, current and density correlations are always intertwined. The latter are sensitive 
to zero-mode contributions from non-zero topological charges. To some extent, it is possible to 
project out these contributions, using 

C lCT2 (i;Pi) = c a airT2 {t ]Pl ) + \a 2 [c^it + l-^ + c^t-l-px)] . (9.3) 

Of course, this projection is only good at small energies in the undesired channel. However, 
we seem to reduce the error in the correlation functions in this way by large factors of about 
20 and more. This may be in part due to the suppression of the zero-mode contributions to 
scalar-density- like correlations, which have particularly high statistical errors (autocorrelation 
times). For the correlations C^ ia2 (t;pi) we obtain sensible fits with the ansatz 

C% 1(J2 (t;p x ) = Y, const -2 cosh[£ CT2 (t - T/2) } . (9.4) 

<T2=± 

The fit interval is chosen so that C" lfT2 (i = 0;pi) does not contribute, except for the correlation 
corresponding to the isoscalar current (//-particle). There the fits are particularly stable with 
respect to the size of the interval size. In this way we always achieve x 2 /d.o.f. < 2. 

9.2 Results 

We now compare results for the n and rj spectrum obtained from standard staggered fermions 
and from our quasi-perfect staggered fermion action. 

For any action based on the truncated perfect free fermion, we expect the dispersion rela- 
tion of the decay energies E{p\) of time slice correlations with fixed spatial momentum p\ to 
approximate the continuum form E(p\) 2 = m 2 + p\- Here, we consider the ^-particle, given by 
the Afjjjl lattice field, and the 7r-particle, given by the M_Z field without disconnected contri- 
butions. The latter choice is motivated by the requirement that masses should be resolved from 
within their error bars, see Section 0. It relies on the absence of disconnected contributions 
to the isovector channel in the continuum. Indeed, for low /3-values, i.e. 7r-masses further away 
from zero, it appears that the masses evaluated from M_Z including disconnected contributions, 
as well as masses from M? , are significantly smaller. Thus we take the worst case of a 7r-mass 
with respect to the continuum massless behavior, and the easiest case for a comparison of the 
different improvement steps. In Figures 12 and 13 we give the energies for the tt and the 77, 
depending on pi for j3 = 3 with the standard staggered action and with the quasi-perfect ac- 
tion. These and the following results are obtained on a 16 x 16 lattice (where we have momenta 
pi = nir/8, n = 0, ...,4). The lines show the continuum dispersion relation with the mass 
derived from p\ = it/ 8. In the quasi-perfect case they fit the data in an excellent way. 



24 




Finally, we show the ^-dependence of the tt- and //-mass in Figures 14 and 15. In the 
Schwinger model, asymptotic scaling for the ry-mass is given by = 2 jn f3. The nearly 
perfect dispersion relation of the quasi-perfect action makes it possible to derive the masses in 
a straightforward manner also from higher momenta, which leads to smaller errors for our data. 

The 77-mass should be equal to the decay mass m ga uge of the plaquette correlation for zero 
spatial momentum. In our case, these are the most precise mass values. They can be extracted 
from a simple cosh- fit. However, for reasons that we do not fully understand, m ga uge does not 



follow the asymptotic scaling prediction as closely as m v , see Figure 16 



10 Conclusions 

Spectral and thermodynamic results for free, truncated perfect staggered fermions have been 



presented before [24], but the improvement did not really reach a satisfactory level there. We 
now pushed that improvement significantly further, mainly thanks to the new blocking scheme, 
which we call partial decimation, but also with the help of a new truncation technique (mixed 
periodic boundary conditions). We now reached a level of excellent improvement, similar to the 
results for truncated perfect Wilson-type fermions. 

We extended the construction of perfect actions to non-degenerate flavors, and we could 
preserve the same level of improvement after truncation also in that case. This is potentially 
important for the study of the decoupling of heavy flavors, or for QCD simulations with realistic 
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Figure 13: The dispersion relation for the n-particle with standard and the quasi-perfect action 
at 13 = 3. 



quark masses. \\ 

At m = 0, which is well described by staggered fermions, the PD scheme is optimally local if 
we just use a 5 function RGT. This simplifies the relation to the continuum n-point functions, 
and numerical RGT steps, which could be performed for interacting theories (in the classical 
limit). 

In d = 2, we added an Abelian gauge field to this fermion along (selected) shortest lattice 
paths plus a "fat link". In the latter, the weight w of the staple term is optimized effectively, 
by minimizing the lowest eigenvalues of the fermion matrix, which led to w = 0.238. 

To simulate the 2-flavor Schwinger model, we applied a variant of the Hybrid Monte Carlo 
algorithm. It uses a simplified action in the identification of possible molecular dynamics steps, 
and the quasi-perfect action in the acceptance decision. In this way, we could avoid an increase 
in the computation effort proportional to the number of couplings. 

As a vertex independent quantity, we measured the Fourier coefficients of the toron distribu- 
tion, and we saw that even the free perfect action (without fat link) helps to move substantially 
closer to the continuum values. 

Next we considered the dispersion relation of the 'V and "77" particle, and we extracted the 
masses m n and at varying (3. We found a tiny pion mass down to f3 ^ 1.5 which confirms 
an excellent scaling behavior of our quasi-perfect action, while m^/S) follows the asymptotic 
scaling very closely. A better scaling is the actual goal of the construction of improved actions, 

7 However, in that case the adequate coupling to gauge fields still has to be worked out. 
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Figure 14: The (J '-dependence of with standard staggered fermions (squares) and quasi- 
perfect fermions, derived from p\ = (crosses) and p\ = tt/8 (diamonds). 



but it was observed before in other models that quasi-perfect actions also tend to improve the 



asymptotic scaling [13, 16]. 



These results, in particular the scaling is even better than the one obtained by using Wilson- 
type fermions with a classically perfect vertex function (parameterized by 123 independent 
couplings) 1 15]. This shows that truncated perfect free staggered fermions, which are suitably 
"gauged by hand", can represent a highly improved, short-ranged action. [] This variant of a 
quasi-perfect improvement program is applicable and promising for QCD. 
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